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The role of Pauli potentials in the scmiclassical simulation of Fermi gases at low tempera- 
tures is investigated. An alternative Pauli potential to the usual bivariate Gaussian form 
by Dorso et aP is proposed. This new Pauli potential allows for a simultaneous good 
reproduction of not only the kinetic energy per particle but also the momentum distri- 
bution and the two-body correlation function. The reproduction of the binding energies 
in finite nuclei in the low and medium mass range is also analyzed. What is found is that 
given a reasonable short-range atractive nuclear interaction one can include correlation 
effects in a suitable chosen density dependent Pauli potential. 
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1. Formalism 



Nuclear many-body simulations are a useful tool to study the relevant properties of 
the nuclear medium in the thermodynamic conditions arising in matter in the after- 
math of a Supernova event or in Neutron Stars. Examples of this are, for instance, 
nuclear pastas 1 ' 2 at densities in the range O.Olpo < P < 0.5po {po = 0.148 /m~ 3 ) and 
temperatures of decens of MeV or in heavy ion collisions. 3 This type of simulations 
based on Monte Carlo or Molecular Dynamics techniques allow for a dynamical de- 
scription of the nuclear medium usually by using an effective interaction hamiltonian 
in a semiclassical treatment. In fermionic systems the genuine antisymmetrization of 
the wave function is considered trough the inclusion of a Pauli potential. Pioneering 
works on this line include those of Wiletsef al. A In this work the hamiltonian used 
to study the low temperature nucleon systems consists of a kinetic energy term and 
a Pauli effective potential (Vp au li)- 
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where 5 TiTj (<5 CTiCTj ) is the Kroneckcr's delta for the nucleon isospin (spin) third- 
component. Pi is the 3- momentum of i-th nucleon and = | — r ^ | (pij — |p, — Pj|) 
the relative distance (momentum) of the i-th and j-th nucleons. 

We will consider, for the sake of comparison, two ways to implement this poten- 
tial. 

i) A Gaussian form introduced by Dorso et al., 5 




Vpaui l {r ljl p lJ ) = V s exp ^~^| - S I > ( 2 ) 

Here po and qo are momentum and length scales related to the excluded phase-space 
volume that is used to mimic fermionic correlations and Vs is the Pauli potential 
strength. All three parameters have been adjusted to reproduce only the kinetic 
energy of a low temperature Fermi gas. 

ii) A new form proposed, based on spatial and momentum-dependent, two-body 
terms of the following form 6 



Vpauu(nj,Ptj) = V q exp(-ry /<7o) + V p exp(-p y -/p ) + V& O n (q t 



(3) 



where % = |p»|/pFj and <d v is a smeared Heaviside-step function, <d v (q) = 
l+cxpt-^g^-i)] anc ^ ®»?(?) — > ®(?) wnen V 1S sufficiently big. The parameters of the 
new Pauli potential V q , V p , Vq and qo,po, r\ will be adjusted to reproduce the kinetic 
energy per particle and both the momentum distribution and two-body correlation 
function of a low-temperature Fermi gas. The first and second terms in the poten- 
tial penalize two particles with the same quantum numbers coming together either 
in space or momentum. This retains the essence of the fermionic wave function 
given by the Slater determinant. The third term forbids any particle from having a 
momentum significantly larger than the Fermi momentum. 

In both cases the potential parameters will depend on the density of the sys- 
tem and this will be crucial in reproducing experimental binding energies when 
simulating low and medium mass nuclei as will be shown later. The values for the 
parameters at saturation density po and reduced temperature r = T/Tp = 0.05 are 
given in Table 1. 



Table 1. Pauli potential parameters. 



Pauli potential 


Potential strength 


90 






(MeV) 


(fm) 


(MeV/c) 


Dorso et al 5 


V s = 207 


1.644 


120 


This work 


V q = 13.517, Vp = 1.260, V e = 3.560 (r; = 30) 


0.66 


49.03 
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2. Results 

The simulations are performed in a NVT system with r = T/Tp, and N fermions 
in a cubic box of volume V = L 3 = N / p. Then, using the Metropolis algorithm the 
system is thermalized until the stage where configurations are sampled in order to 
calculate the statistical averages for the magnitudes discussed below. 




Fig. 1. Energy per particle of a Fermi gas simulated with N = 1000 particles at T = 0.05Tp as 
a function of density. 

In Fig.[T]the red line shows the kinetic energy per particle for a Fermi gas system 
with TV = 1000 particles at r = 0.05 as a function of density calculated using the 
new form of the Pauli potential Eq.Q. Also plot with a black line is the exact 
result. We can see that there is a good reproduction of the kinetic energy. 




Fig. 2. Momentum distribution function and two-body correlation function of a Fermi gas. See 
text for details. 

In Fig. [2] we can see on the left side the momentum distribution function f(q, r) 
and, on the right side, the two-body correlation function g{z,r) with z — pft. The 
blue curves correspond to the Dorso potential Eq.([2]) and the red curves to the 
new Pauli potential proposed Eq.Q. Again the black line shows the exact result. 
We can see that a simultaneous good reproduction of both magnitudes is achieved 
with the alternative new potential but not with the Dorso version. Particurlarly the 
"Fermi hole" fails to be reproduced at small distances with the Dorso potential. This 
should be emphasized since these models are used in nuclear many-body simulations 
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Fig. 3. Application to nuclei, (a) Binding energy, (b) Kinetic and potential contributions to the 
binding energy. 



as in nuclear pastas as, for instance, in the work by Maruyama et al. 2 The velocity 
distribution, 6 not shown here, however peaks at lower values than the momentum 
distribution due to the fact that canonical and kinematical momentum are not 
the same quantities. 7 This is a genuine feature in this treatment with momentum 
dependent Pauli potentials in a hamiltonian formalism. 

We now show finite nuclei simulation 8 results calculated with a simplified square- 
well nuclear potential with V we ii = —3 MeV of width 2 fm and a core with V core = 10 
MeV and width 1 fm. Coulomb interaction is also included. In Fig. 02a) binding 
energy per particle for a low to medium mass set of spin saturated symmetric nuclei 
oi A nucleons. As can be seen in Fig. [Hb) kinetic (dashed line) and potential(dotted 
line) energy balance to obtain the total binding energy (solid line) per particle. The 
density dependence of the parameters of the Pauli potential is crucial to provide 
enough positive contribution to the linearly A-growing negative potential energy 8 
and reproduce the experimental binding energy curve. 
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